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ABSTRACT 

Coronal jets are collimated, dynamic events that occur over a broad range of spatial scales in 
the solar corona. In the open magnetic field of coronal holes, jets form quasi-radial spires that 
can extend far out into the heliosphere, while in closed-field regions the jet outflows are confined 
to the corona. We explore the application of the embedded-bipole model to jets occurring in 
closed coronal loops. In this model, magnetic free energy is injected slowly by footpoint motions 
that introduce twist within the closed dome of the jet source region, and is released rapidly by 
the onset of an ideal kink-like instability. Two length scales characterize the system: the width 
(N) of the jet source region and the footpoint separation (L) of the coronal loop that envelops the 
jet source. We find that the jet characteristics are highly sensitive to the ratio L/N, in both the 
conditions for initiation and the subsequent dynamics. The longest-lasting and most energetic 
jets occur along long coronal loops with large L/N ratios, and share many features of open-field 
jets, while smaller L/N ratios produce shorter-duration, less energetic jets that are affected by 
reflections from the far-loop footpoint. We quantify the transition between these behaviours 
and show that our model replicates key qualitative and quantitative aspects of both quiet-Sun 
and active-region loop jets. We also find that the reconnection between the closed dome and 
surrounding coronal loop is very extensive: the cumulative reconnected flux at least matches the 
total flux beneath the dome for small L/N, and is more than double that value for large L/N. 

Subject headings: Sun: corona; Sun: magnetic fields; Sun: jets; magnetic reconnection 


1. Introduction 

Observations of the Sun’s outer atmosphere - 
the chromosphere, transition region, and corona 
- reveal the ubiquitous occurrence of jetting phe¬ 
nomena across a wide range of spatial and tempo¬ 
ral scales. Transient, impulsive, collimated flows 
of plasma are observed both as bright, emitting 
features at high (coronal) temperatures and as 
dark, absorbing features at low (chromospheric) 
temperatures; generally, these events are referred 


to as jets and surges, respectively (e.g. Canfield 


et al.|1996 ). The wavelengths involved range from 
the optical, through the UV and EUV, to X-rays. 
Ha surges were identified and studied first, from 
ground observatories (e.g. Roy||1973 ). The advent 
of telescopes placed in space, above Earth’s at¬ 
mosphere, was a prerequisite for the detection of 


UV and EUV jets from Skylab (Schmahl 1981) 


and from sounding rockets (Brueckner & Bartoe 


1983), and, later still. X-ray jets from Yohkoh 


(Shibata et al. 1992, 1994 Shimojo et al. 1996). 


Cool surges and hot jets sometimes are observed 
together, in close association in space and time 
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( Canfield et aL|p!9% ). Subsequent improvements 
in the spatial resolution and temporal cadence of 
space-borne instruments have enabled ever more 
detailed studies of jets and surges from the Solar 


and Heliospheric Observatory {SOHO) (e.g. Wang 
et al.|p!9M ), Transition Region and Coronal Ex¬ 


plorer {TRACE) (e.g. |Chae et al.| 1999[ ), Solar- 


Terrestrial Relations Observatory {STEREO) (e.g. 


Patsourakos et al.j2008|), Hinode (e.g.|Cirtain et al. 

to 

o 

o 

Savcheva et al.||2007t |Nishizuka et al.||2008t 

Torok et al. 2009| Moore et al.||2010 

2013| |Lm| 

et al. 2011), Solar Dynamics Observatory {SDO) 

(e-g- 

Guo et al. 

2013 Lee et al. 2013 

Schmieder 

et al. 2013 Zheng et al. 

2013 

Shen et al.| 2011|), 


and Interface Region Imaging Spectrograph {IRIS) 
(e.g. Tian et a'r]|2Q14 Cheung et al.||2Ql^ . 


Typically, an EUV or X-ray jet begins with a 
rapid brightening low in the solar atmosphere, in¬ 
dicating an impulsive increase in the plasma tem¬ 
perature or density, or both. This is followed 
by bulk outflows of material from the bright re¬ 
gion, usually at supersonic speeds and highly colli¬ 
mated in direction. These properties suggest that 
the plasma flow is guided along the local mag¬ 
netic field - tracing loops in magnetically closed 
active regions and quasi-radial spires in magnet¬ 
ically open coronal holes - and is subjected to 
nonthermal, presumably magnetic, forces. The 
bright emissions fade away gradually as the plasma 
cools and expands on its passage through the at¬ 
mosphere and as the energy source driving the 
jet is depleted. In many jets, the outflows ex¬ 
hibit a distinctly helical structure in the loop or 
spire. Within this subclass, many events also dis¬ 
play translational bodily motions of the jet loop 
or spire across the plane of the sky. With or with¬ 
out these helical and translational motions, a mi¬ 
nority of jets are observed to recur from a single 
structure, yielding multiple episodes of brighten¬ 
ing, jetting, and dimming that are separated by 
intervals of quiet. 


All of the above properties have been es¬ 
tablished by numerous reported observations of 
coronal-hole jets, which frequently are referred 
to as polar jets. These characteristics also have 
been observed, albeit for a much smaller event 
sample, in jets occurring in closed coronal loops 


(e-g- 

Torok et al.|2009 Yang et al.|2012||Guo et al. 

2013 

Lee et al.|2013t Schmieder et al.|2013| Zheng 

et al. 

2013 Gheung et al. 2015). The magnetic 



Fig. 1.— A closed-loop jet observed by SDO, de¬ 
scribed in detail by Cheung et al.| ( 2Q15| ). HMI 
line-of-sight magnetograms (grayscale) overlaid 
with AIA 94 A and 131 A channel images (blue) 
show the (a) launch and (b) propagation to the 
far footpoint of a helical coronal-loop jet. An ani¬ 
mation of this figure is available online. 


structures hosting these jets range in size from 
long loops rooted in weak-field areas adjacent to 
coronal holes, to short loops rooted in strong-held 
areas within young active regions. One example 
of a closed-loop jet, observed with SDOA Atmo¬ 
spheric Imaging Assembly (AIA) and described in 
detail by Cheung et al. I dM^ , is shown in Figure 
Our objective in this paper is to simulate jets 
in a variety of such closed coronal loops by em¬ 
ploying a well-tested model, described below, that 


has been applied extensively to polar jets ( 

Pariat 

et al. 

120091 |2010| IDalmasse et al.||2012[ 

Pariat 

et al. 

2015). 


The source region of many solar jets has the 


morphological appearance of a sea anemone (Shi- 


bata et al. 1994), with bright tendril-like curved 
loops emanating radially from a central locus. It 
is well established from observations that such 
anemone regions form as a result of hux emer¬ 


gence into coronal holes (e.g. Torok et al. 2009 Liu 
eFaLllMI ). The photospheric magnetic held be- 
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neath the loops typically has a concentrated patch 

Shibata||1995||1996t |Miyagoshi &; Yokoyama 

2003 

of vertical field of one (minority) polarity embed- 

2004| Archontis et al.||2005| |Galsgaard et al. 

2005 

ded within a large-scale area of generally weaker 

Moreno-Insertis et al.|2008||Gontikakis et al. 

2009 

vertical field of the other (majority) polarity (e.g. 

Torok et al.|2009t Archontis et al. 

20101 Jiang et al. 


Cheung et al. 2015). Magnetic flux tubes that 


originate within the minority-polarity patch close 
back to the Sun’s surface locally, on the far side 
of the polarity inversion line encircling the patch. 
The manner in which the two polarity regions are 
connected depends on the free energy stored in 
this closed field. The minimum-energy potential 
magnetic field above such a distribution of photo- 
spheric flux has a dome-shaped structure with a 
magnetic null point near its top. Flux tubes in¬ 
side of the dome form closed, unsheared loops like 
those observed in the anemone regions; flux tubes 
outside of the dome follow the large-scale back¬ 
ground field associated with the majority polarity, 
closing back to the Sun farther away (in the case 
of closed-loop jets, as in the observational papers 
cited earlier and as shown in Fig.[^ or opening out 
into the remote heliosphere (in the case of coronal- 
hole jets). Two of our model setups exhibiting a 
potential null dome embedded within closed coro¬ 
nal loops are shown in Figure The separation of 
the magnetic configuration into distinct flux sys¬ 
tems, one closed locally and the other open (or 
closed remotely, as in Fig. |^, readily allows rel¬ 
ative displacements of field lines to occur across 
the null point. Such displacements generate strong 
electric currents and, eventually, initiate magnetic 


reconnection and associated activity (Lau & Finn 
!^tiochos|p!9^ ). 


1990 


A prototypical model for solar jets, account¬ 
ing for their impulsiveness and their helical mo¬ 
tions by appealing to magnetic reconnection be¬ 
tween a closed volume of twisted magnetic flux 
and the ambient open untwisted field, was put 
forth by Shibata & Uchida (1986). The reconnec¬ 


tion transfers magnetic twist from closed to open 
field lines, and the twist residing on the newly re¬ 
connected open field lines then propagates away 
from the interaction region at the Alfven speed to 
form an “untwisting” jet. This basic picture has 
been the basis for numerous simulations of jets 
driven by flux emergence through the solar pho¬ 
tosphere, in which a twisted flux rope rises buoy¬ 
antly into a pre-existing coronal magnetic field, 
initiates reconnection, and launches a jet when 
the field orientations are favorable (Yokoyama & 


2Q12| Archontis fc Hood||2Q13 N oreno-Insertis fc 


Galsgaard||2Q13| [Takasao et al.|[2QT3l [Fang et ah 

2Q14| ). This process also can form the anemone 


structure itself, as the observed end state of the 
evolution following closed-loop ( Torok et al.|2QQ9 ) 
and coronal-hole ( Liu et a'r]|2Qll ) jets. 

It is by no means clear that significant flux 
emergence always precedes solar jets, however. 
The fundamental susceptibility of the null-point 
configuration to reconnection can be exploited in 


other ways, as noted by Shibata & Uchida (1986). 


In studies of polar jets in open fields, Pariat and 
collaborators ( Pariat et al.||2QQ9 2010 |2015l |Dal- 


masse et al. 2012) initiated reconnection and jet 
onset by slowly twisting the closed flux beneath 
the anemone dome with applied footpoint mo¬ 
tions. When the amount of twist reaches a critical 
threshold, which depends upon the inclination of 
the ambient field to the surface, the closed flux 
convulses due to onset of an ideal kink-like insta¬ 
bility ( Rachmeler et al.||2010 ). The kink itself re¬ 
leases very little energy, but it initiates impulsive 
reconnection between the twisted closed flux inside 
and the untwisted open flux outside of the dome. 
This reconnection rapidly and efficiently releases 
a large fraction of the stored magnetic free energy, 
launching a very strong helical jet. In addition, 
as the reconnected open field “untwists” it carries 
the reconnection site around the dome, generat¬ 
ing a translational bodily motion of the resulting 
jet spire. The reconnection and the jet flows sub¬ 
side gradually as the stored magnetic energy is 
depleted. If, however, the footpoint motions per¬ 
sist, then the free energy builds up again until an¬ 
other cycle of reconnection, jetting, and relaxation 
occurs, yielding recurrent events. All of these fea¬ 
tures are consistent with observed EUV and X-ray 
jets. 

In this paper, we explore the ramifications for 
the embedded-bipole model of jets when the am¬ 
bient magnetic structure is a closed coronal loop, 
rather than an open coronal hole. As may be an¬ 
ticipated, the behavior is quite similar to the polar 
case for long loops; for short loops, on the other 
hand, there are both quantitative and qualitative 
differences in the dynamics. We describe our re- 
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suits in the rest of the paper, which is organized 
as follows. In ^ we describe the physical and nu¬ 
merical model, ^compares and contrasts two rep¬ 
resentative cases, a long coronal loop and a short 
one, while il describes the broader results of a 
parametric survey over loop length. We relate our 
simulations qualitatively and quantitatively to jet 
observations in il A summary of the highlights of 
our results is given in il along with a discussion 
of their implications for future research. 

2. Model 

We consider the simplest possible magnetic con¬ 
figuration for a jet-generating anemone region, in 
which a small-scale patch of relatively strong ver¬ 
tical field of one polarity is embedded in a large- 
scale region of weaker vertical field of the opposite 
polarity. This large-scale field, in turn, belongs 
to a bipolar flux distribution that forms a long 
coronal loop. Two particular realizations of this 
configuration are shown in Figure The topol¬ 
ogy is of the fan/spine type, with a magnetic null 
point atop an approximately hemispherical dome 
of fan field lines. Emanating from the null are 
inner and outer spine lines that root in the par¬ 
asitic polarity and the far polarity of the back¬ 
ground field, respectively. Two intrinsic length 
scales of this system are the diameter {N) of the 
separatrix dome and the separation (L) of the two 
spine foot points. Both of these lengths are easy 
to measure numerically, and they can be approx¬ 
imated from observational data as the width of 
the brightened anemone region and the footpoint 
separation of the large-scale coronal loop enclosing 
the anemone, respectively. The aspect ratio L/N 
of the configuration quantifies the relative sizes of 
the coronal loop and the enclosed separatrix dome, 
ranging from nearly unity when the background 
loop is almost as compact as the dome itself to 
very large values when the background loop is far 
larger in scale. Open-field jets correspond to the 
limit L/N oo, where the “coronal loop” extends 
from the Sun out to the remote heliosphere. 

The scenarios constructed and simulated nu¬ 
merically for this paper were all performed in non- 
dimensional units. This is convenient, given the 
broad ranges of scales of jets and coronal loops 
on the Sun. It is allowed due to the homogene¬ 
ity of the magnetohydrodynamics (MHD) model. 



- 4.0 2.5 9.0 15.5 22.0 

Bx 


Fig. 2.— Initial magnetic field in two config¬ 
urations with aspect ratios L/N = 2.40 (top) 
and L/N = 1.46 (bottom). The bottom planes 
are color-shaded according to the sign (+, —) and 
strength of the field component normal to the sur¬ 
face (Bx). Selected magnetic field lines outline 
the fan separatrix surface and the inner and outer 
spine lines emanating from the null point. The di¬ 
ameter of the separatrix dome (N) and the sepa¬ 
ration of the spine-line footpoints (L) characterize 
the intrinsic spatial scales of the configuration. 
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which permits fundamental scales for length, time, 
and mass to be factored out of the set of variables 
and equations. All of the dimensionless numbers 
governing the physics of the system - the acoustic 
and Alfven Mach numbers, plasma beta, Reynolds 
number, Lundquist number, etc. - are preserved 
under such rescalings of the fundamental length, 
time, and/or mass. Later in ^5.2, we will scale our 
non-dimensional results by assuming some typi¬ 
cal values for solar parameters, thereby deriving 
predicted properties of the simulated jets for com¬ 
parison with observations. 


2.1. Initial Conditions 

The magnetic configurations illustrated in Fig¬ 
ure are minimum-energy, electric current-free 
(“potential”) fields. To construct them, we super¬ 
posed two magnetic dipoles: one horizontally ori¬ 
ented to generate the large-scale background field; 
the other vertically oriented with a more compact 
flux distribution at the coronal base to generate 
the parasitic polarity. By changing the strength 
and position of the vertical dipole, we varied the 
location and relative size of the separatrix dome 
compared with the coronal loop that encloses the 
dome. 

Specifically, the initial magnetic field is given 
by 

B = V X (A/i + Ay), (1) 

where and are the vector potentials for the 
horizontal and vertical dipoles respectively. We 
orient the x direction vertically in our simulation 
domain, with x = 0 at the coronal base (the “pho¬ 
tosphere”), so the vector potentials can be written 

^ ^ Bhdj _ {z - Zh)±- {x - Xh)z _ 

2 [{x - Xh)'^ + {y-yh)'^ + (z - 

^ Bydl -{z - Zy)y + {y - yy)z 

2 [{x - Xy)"^ + {y - yy)"^ + {z - Zy)‘^f^‘^ 

(2) 

We placed the background horizontal dipole below 
the photosphere at {xh^yh^^h) = (—10.0,0.0,0.0), 
so its depth dh = 10.0, with signed strength 
Bh = +8.0. This produces a peak field strength 
\'^h\max = 4.0 at the photosphere above the hor¬ 
izontal dipole, {x^y^z) = (0.0, 0.0, 0.0). The ver¬ 
tical dipole was placed below the photosphere at 
{xy,yy,Zy) = (-1.7, ^^,0), SO Its depth dy = 1.7, 


with signed strength By = +25.0. We varied the 
y position of the vertical dipole from yy = —4 to 
yy = —10. Its orientation is such that it opposes 
the direction of the background dipole field at its y 
position, generating the parasitic polarity and cre¬ 
ating the 3D null-point topology. Figure shows 
magnetic field lines in the z = 0 plane and depicts 
the positions of the two magnetic dipoles relative 
to the photosphere for the case L/N = 2.40. The 
peak vertical magnetic field above the parasitic po¬ 
larity varied little {B^^^ 21-22) over the range 

of yy values considered. The parameters of our 
various calculations are given in Table These 
include yy, N, L, L/N, and the positive magnetic 
flux under the separatrix dome, "^dome- 

The peak magnetic pressure associated with our 
large-scale background field is |B/i|^^^/87r 0.64. 

We chose P = 0.01 for the dimensionless thermal 
pressure, yielding a minimum plasma beta (ratio 
of thermal to magnetic pressure) of jSh ~ 1.5 x 10“^ 
with respect to the background field. At the 
center of the parasitic polarity, where the field 
strength is higher, the plasma beta is lower still, 
Py Qx 10“^. The spatially varying field strength 
means that near the photosphere (3 generally is 
well below unity and the magnetic field strongly 
dominates the dynamics, whereas high in our 
gravity-free corona, where the field is weaker but 
the thermal pressure remains uniform, (3 reaches 
and then slightly exceeds unity. At and near the 
null point, in contrast, (3 becomes very large (infi¬ 
nite, in principle, at the null) and the plasma dom¬ 
inates the dynamics, as it should. The majority 
of the jets in our parametric survey remain con¬ 
fined within the \ow-P region; however, the exper¬ 
iments with > 8.5 produce jets that impinge 
on the moderate-/3 high corona. In these jets, the 
outflows reaching the apex of the loop bend the 
held lines rather sharply in the moderate-beta re¬ 
gion. This bending of the held does not impede 
the progress of the jet, however, and the moder¬ 
ate (3 in this region has no effect on the jet trigger 
and energy release, which occur well down in the 
low-/d portion of the corona. 

We can set the mass density independently of 
the thermal pressure, so we chose an initial uni¬ 
form value p = 1. The selection R = 0.01 for the 
ideal gas constant then sets the units of tempera¬ 
ture, with an initial uniform value T = 1. These 
additional choices result in a uniform sound speed 
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4.0 

4.5 

5.0 

5.5 

6.0 

6.5 

7.0 

7.5 

8.0 

8.5 

9.0 

10.0 

N 

5.86 

5.82 

5.82 

5.84 

5.88 

5.96 

6.08 

6.18 

6.34 

6.48 

6.69 

7.04 

L 

6.02 

7.32 

8.50 

9.68 

10.82 

11.96 

13.06 

14.16 

15.20 

16.26 

17.31 

19.24 

L/N 

1.03 

1.26 

1.46 

1.66 

1.84 

2.01 

2.15 

2.29 

2.40 

2.51 

2.59 

2.73 

^ dome 

52.1 

51.3 

50.9 

51.0 

51.5 

52.2 

53.2 

54.3 

55.5 

56.9 

58.3 

61.2 

rptot 

-^inj 

115.5 

118.1 

120.3 

122.0 

123.9 

124.6 

126.4 

124.6 

131.9 

132.7 

133.4 

134.5 

ttrig 

1180.0 

1010.0 

840.0 

790.0 

720.0 

655.0 

535.0 

560.0 

570.0 

650.0 

790.0 

- 

tjet 

80.0 

130.0 

100.0 

130.0 

120.0 

165.0 

260.0 

300.0 

315.0 

330.0 

330.0 

- 


Table 1: Simulation parameters: Vertical dipole position (l^^l), separatrix dome diameter (V), length of 
outer coronal loop (L), aspect ratio (L/V), positive magnetic flux under the separatrix dome {"^dome)^ total 
injected energy and measured jet trigger times and durations {Urig and tjet, respectively). 



Fig. 3.— Magnetic field lines (solid curves) in the 
z = 0 plane show the configuration with L/N = 
2.40 constructed using two sub-photospheric mag¬ 
netic dipoles. The dashed line shows the position 
of the photosphere, which is color-shaded in the 
top panel of Figure 2 . 


Vs ~ 0.13 throughout the domain, for a ratio of 
specific heats 7 = 5/3. The same parameter val¬ 


ues were assumed by Pariat et al. (2009 2010 


2015). The maximum Alfven speed is Va ~ 5.9 
at the center of the parasitic polarity. At the 
mean value across our sample of the aspect ratio, 
L/V ^ 1 . 88 , the diameter of the closed separatrix 
dome is V 5.9. Consequently, one dimensionless 
time unit in our simulations equals one character¬ 
istic Alfven time defined by these values for Va and 
N. 


2.2. Temporal Evolution 


We use the Adaptively Refined Magnetohydro¬ 


dynamics Solver (ARMS; DeVore & Antiochos 
2008) to solve the ideal MHD equations in the form 

| + v^(^) = o, 

(V X B) X B 


(3) 


dt 


V • (pvv) + VP — 


/io 


dU 

dt 


^-h V • (Pv) + PV • V = 0, 


— - V X (v X B) = 0 . 


= 0 , 

(4) 

(5) 

(6) 


Here t is the time, p the mass density, P = pRT 
the thermal pressure, U = P /(7 — 1 ) the inter¬ 
nal energy density, /io = 47 r the magnetic per¬ 
meability, and B and v the 3D magnetic and ve¬ 
locity fields. Reconnection occurs in our simu¬ 
lations through numerical diffusion terms in the 


Flux-Corrected Transport scheme (DeVore 1991) 


employed by ARMS. The lack of an explicit re¬ 
sistivity allows the greatest amount of free mag¬ 
netic energy for the given resolution to be built 
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up in our calculations before rapid reconnection 
is initiated and the jet occurs. The FCT solution 
algorithm adds explicit numerical diffusion to the 
equations, and a corresponding amount of explicit 
anti-diffusion that minimizes the residual trunca¬ 
tion error in smooth regions of the flow. These 
error-cancelling anti-diffusion fluxes are limited 
(‘corrected’), however, so that artificial numerical 
ripples are suppressed at shock fronts, shear lay¬ 
ers, current sheets, and other discontinuities, as 
described in detail by DeVore (1991). The effect 
is to introduce just enough dissipation to keep the 
solution well-behaved when ideal motions, such as 
the onset of kink instability in the case of our jet 
simulations, induce sudden, strong changes in the 
magnetic field and the plasma variables. 

We energize the system in the same manner 
Pariat et al. (2009). That is, we impose a 


as 


slow, subsonic and sub-Alfvenic, rotation of the 
parasitic polarity that introduces magnetic shear 
across the nearly circular polarity inversion line. 
Specifically, the driving profile is 


v_L = vof{t)g{Bx)± X VBx, 


m ^ 2 


1 — COS 27r 


tt 


g{Bx) = ks tanh (^b 


Bx - Bi 
Br — Bi 


(7) 


The above expressions are used for t G [0, ttwist] 
and Bx G outside of those intervals, we 

set = 0. The flow follows the contours of Bx 
within the parasitic polarity patch and is incom¬ 
pressible, i.e. divergence-free, so that it preserves 
the vertical component of the magnetic field at 
the photosphere throughout the evolution. The 
flow speed vanishes exactly as Bx Bi^ and it 
becomes small as Bx B^ because \/Bx becomes 
small at the center of the parasitic polarity patch. 
We chose Bwist = 1000, Bi = 0.6, ks = 5.0, and 
vq = 6.84 X 10“^; B^ was varied over the narrow 
range [20,21], depending upon the peak strength 
of the dipole, to rotate the dipole as close to its 
central axis as possible. For these choices, between 
0.7 and 1.1 maximum turns of twist, M, are in¬ 
jected into the field beneath the dome by the time 
of jet initiation, as shown below. The peak driving 
speed on the boundary is |vx| ~ 0.016, which is 
about 12% of the sound speed and 0.3% of the lo¬ 
cal Alfven speed. Thus, the magnetic field evolves 


quasi-statically and remains approximately force- 
free throughout the low-/3 portion of the corona. 

The adaptive grid employed by ARMS is con¬ 
structed from a basis set of root blocks, which can 
be subdivided to attain higher grid refinements in 
a pre-defined way and/or adaptively as the solu¬ 
tion requires ( MacNeice et al.||2000 ). The root 
blocks in these simulations had a fixed spatial ex¬ 
tent of 17 X 17 X 17. For values of \yy\ in the range 
[4, 8], a domain size of [0, 34] x [—17,17] x [—8.5, 8.5] 
(2x2x1 blocks) was sufficient to avoid any signifi¬ 
cant influence of the boundaries during the jet evo¬ 
lution. For larger values of |^^|, a larger domain 
was necessary. To maintain the same grid spacing 
in these calculations, the number of root blocks 
was increased to 2 x 3 x 2, giving a domain size of 
[0,34] X [—25.5,25.5] x [—17,17]. In each calcula¬ 
tion, we required a minimum grid refinement of 4 
and a maximum of 6 from the initial 32 x 32 x 16 
grid (on the smaller domain). A volume in each 
simulation that encompasses the footprint of the 
separatrix dome, with a height of ^ 0.5, is fixed 
at the maximum refinement level to resolve the 
boundary driving motions as finely as possible. 
The grid outside this volume adapts according to 
whether strong gradients, beyond a fixed dimen¬ 
sionless value, develop in the magnetic field (for 
details, see Karpen et al. 2012). This adaptive 
refinement better resolves the electric current lay¬ 
ers and any shocks that develop in the domain. 
For a typical calculation on the smaller domain 
this results in approximately 1.2 x 10^ grid cells, 
compared with around 5.4 x 10^ for the equiv¬ 
alent grid uniformly refined to the same maxi¬ 
mum resolution everywhere. The larger domain 
required only slightly more grid cells, because the 
outer regions where the domain was extended were 
refined almost solely to the minimum level of 4. 
All boundaries of the box were closed (zero fluxes 
of mass, momentum, energy, and magnetic flux 
passed through) and line-tied (v = 0 except where 
nonzero was imposed according to Equation]^. 


3. Results I: Long and Short of Loop Jets 

We present our results by first focusing on the 
differences between the jets produced in the two 
configurations shown in Figure a short coro¬ 
nal loop with small aspect ratio, L/N = 1.46, and 
a long loop with large aspect ratio, L/N = 2.40. 
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The latter value is large enough to exhibit dynam¬ 
ics that are markedly different from the former, 
while also resembling rather closely those of the 


open-field jets simulated by Pariat et ah (2009 


2010 2015). These two cases fall near the ex¬ 


tremes of our parameter range, which will be ex¬ 
plored more fully below in Q The variation in 
L/N is determined principally by the coronal loop 
length, L G [6.0,19.2], whereas the size of the sep- 
aratrix dome is roughly fixed, N G [5.8, 7.0]; see 
Tabled 

There are two aspects to the dependence upon 
L/N. First, the explosive reconnection that gives 
rise to a jet is of finite duration, as is the time 
required for disturbances launched by the jet to 
reach the far footpoints of the coronal loop and 
refiect back into the jet source region. For long 
loops, the duration is less than the travel time, 
and the jet generation process essentially is unaf¬ 
fected by the fact that the enveloping structure 
closes back to the Sun. For short loops, on the 
other hand, the duration can be longer than the 
travel time, and refiected disturbances can impact 
the jet source region before the generation process 
is complete. Our two examples illustrate this dis¬ 
tinction. Second, the position of the null point 
on the separatrix dome is sensitive to the incli¬ 
nation of the background horizontal dipole field 
above the position of the embedded polarity. For 
long loops, the orientation is nearly vertical, the 
null point is positioned near the top of the dome, 
and the entire dome is nearly axisymmetric. For 
short loops, in contrast, the orientation is far from 
the vertical, the null point is positioned well over 
on the side of the dome nearest the far loop foot- 
point, and the dome is strongly asymmetric. In 
the latter case, it is much easier to displace the 
inner and outer spine lines along the fan surface, 
thereby distorting the potential null point into a 


current patch where reconnection can occur (An- 


tiochos 


1996 Pontin et al. 2007). In the former 


case, the reconnection is impeded very effectively 
until onset of an ideal kink-like instability strongly 
breaks the near-axisymmetry of the dome. A sur¬ 
vey of simulated open-field jets with varying tilt 


angle of the uniform background field (Pariat et al. 


2015) confirms this expected range of behaviors. 


We anticipated that our more asymmetric, short- 
loop configuration would form a current layer more 
readily and that reconnection would play more of 


a role at all stages of the evolution than for our 
more symmetric, long-loop configuration. As de¬ 
tailed below, this is just what we observe. 


(2009 2010 2015) 


In all of the open-jet calculations of [Pariat et al.| 
the system evolved through 
Our new closed-field simula- 
During the ini- 


three main phases, 
tions also exhibit this progression 
tial, energy-storage phase, the twist slowly builds 
in the field beneath the separatrix dome and mag¬ 
netic free energy accumulates in the structure. 
There is little to no energy release due to any slow 
reconnection that occurs at the strengthening null- 
point current patch. Eventually, a critical thresh¬ 
old for magnetic twist or free energy - governed by 
the ideal kink-like instability for the axisymmetric 
case, at least - is attained, the separatrix dome 
convulses, and the impulsive energy-release phase 
begins. Rapid spine-fan reconnection ([Priest & 


Pontin[[2QQ9 ) is initiated across the separatrix, re¬ 


leasing much of this twist and free energy onto 
external, shear-free field. The untwisting of the 
newly reconnected field lines causes the reconnec¬ 
tion site to process around the separatrix dome, 
generating nonlinear torsional Alfven waves that 


comprise a helical jet (Patsourakos et al. 2008 


Pariat et al. 2009). Subsequently, these waves 


propagate away from the separatrix dome during 
a concluding relaxation phase, transporting signif¬ 
icant magnetic energy and helicity away from the 
jet source region and causing the jet to subside as 
the field relaxes toward a lower energy state. 


The following subsections qualitatively com¬ 
pare and contrast the three phases of jet evolu¬ 
tion for our two selected cases. A final subsection 
analyzes the rate and location of the reconnection 
that occurs in these events before a quantitative 
discussion of the full parameter study in Q 


3.1. Energy Storage Phase 

The evolution of the system as twist is intro¬ 
duced by the boundary driving motions prior to 
the jet is shown in Figure^ For the long-loop 
case {L/N = 2.40, Fig. |^-c), the closed fiux 
within the separatrix dome (yellow field lines) ex¬ 
pands upward as the magnetic pressure increases 
beneath the dome. A current layer of small spa¬ 
tial extent gradually forms about the null as the 
dome expands into the surrounding loop fiux (Fig. 

-c). The reconnection associated with this cur¬ 
rent layer is small, and the outfiows are weak, so 
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Fig. 4.— Energy buildup phase for L/N = 2.40 (top) and 1.46 (bottom). Yellow and purple magnetic field 
lines are traced from fixed, line-tied footpoints that at t = 0 reside inside and outside of the separatrix, 
respectively. The horizontal-plane shading shows as in Figure The vertical-plane contours show the 
current-density magnitude |J|, whose scale is saturated for clarity. Animations of both configurations are 
available online. 


that very little of the flux beneath the dome re¬ 
connects prior to jet onset. In contrast, for the 
short-loop case {L/N = 1.46, Fig. [^-f), a much 
more extended current layer quickly forms around 
the null and the nearby separatrix surface as the 
dome expands. The resulting reconnection links 
flux previously closed beneath the dome to the far 
coronal loop footpoint (Fig. [^, yellow field lines), 
thereby transferring part of the injected twist onto 
the enclosing loop. Indeed, in §4.2| we will show 
that almost all of the flux beneath the dome is re¬ 
connected during this phase prior to onset of the 
short-loop jet. 

3.2. Energy Release Phase 

In both cases, the initiation of fast energy re¬ 
lease appears to be driven by the onset of a kink¬ 
like instability. Figure]^ shows the field line (yel¬ 
low) with the greatest number of turns beneath 
the dome prior to the initiation of the impulsive 
jet in each case {Mfrig ~ 0.8 and 1.1 turns 
for L/N = 2.40 and 1.46, respectively). Also 
shown are field lines (blue) traced from evenly 
spaced footpoints around the contour of from 
which the yellow field line begins. Because our 
boundary driving follows the contours of Bx^ al¬ 


beit at non-constant speed, under ideal evolution 
the blue and yellow field lines should have approx¬ 
imately the same number of turns about the inner 
spine of the null in each simulation. For the long 
loop with L/N = 2.40, this is indeed the case, 
showing that the weak reconnection near the null 
has had little effect on the most sheared magnetic 
flux. The approximate cylindrical symmetry of the 
most strongly sheared field is maintained up to the 
time of initiation of the jet. The sudden break¬ 
ing of the symmetry and onset of impulsive recon¬ 
nection closely resembles the evolution of previous 
open-field calculations with perfect initial symme¬ 


try ( Pariat et al.|2QQ9 Rachmeler et al.|2QlQ ), and 
is consistent with onset of a kink-like instability. 
The critical number of turns at onset, Mtrig ~ 0.8, 
is very close to that reported by Pariat et al. (2010) 


when the background open field is tilted from the 
vertical by 10°. Our long-loop configuration has 
an effective tilt angle of about 9° (see ^4.1 and Fig. 


II 


for more details). 

For the short loop with L/N = 1.46, in contrast 
to the long loop, the early reconnection reaches 
field lines that are close to the PIL, along most 
of its extent. This acts to reduce the shear of the 
field lines that straddle the PIL so that some field 
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Fig. 5.— Injected twist just prior to jet initiation. 
The magnetic field line with the highest number 
of turns about the inner spine is colored yellow. 
Blue field lines are traced from footpoints evenly 
spaced along the contour of (horizontal-plane 
shading) where the yellow field line starts. Silver 
field lines show the field structure near the null 
region. 


lines on the same contour of are more sheared 
than others. Along a given contour, the unrecon¬ 
nected field (having the greatest number of turns) 
folds underneath these less-twisted field lines (Fig. 

yellow field line), thereby contorting the shape 
of the separatrix. Thus, the early reconnection in 
the short-loop case reduces the amount of flux be¬ 
neath the dome that is most strongly twisted prior 
to initiation of the instability. It also enhances the 
cylindrical asymmetry of the configuration. Nev¬ 
ertheless, the maximum number of turns attained 
{Mfrig ~ IT) is well above that achieved (~ 0.8) 
prior to the long-loop jet, while still well below 
that of the perfectly symmetric open case (~ 1.4). 
This is also true of the average number of turns 


( p.lj and Fig. 10). Evidently, the critical twist for 
onset of the kink-like instability driving these im¬ 
pulsive jets is not a simple, monotone function of 
the effective tilt angle of the configuration. 

The jet produced in our long loop {L/N = 2.40; 
Fig.|^-c) is qualitatively similar to those observed 
in the open configurations of Pariat et al. ( |2009 


The jet starts at ttrig ~ 570 and 
lasts for tjet ~ 315 Alfven times. The jet is initi¬ 
ated when the instability onset forcibly reconfig¬ 
ures the sheared configuration, driving the twisted 
field into the underside of the separatrix dome and 
generating an extended helical current layer across 
which magnetic flux reconnects rapidly. The fast 
reconnection-driven plasma outflows that follow 
this transition are shown as isosurfaces of veloc¬ 
ity magnitude in Fig.the isosurfaces are colored 
according to the local value of Vz^ showing the ro¬ 
tational component of the flows as they progress 
along the loop. Following initiation (Fig. [^), as 
the reconnection continues and the helical current 
layer forms, the outflow speed increases and takes 
on a more helical shape (Fig. |^). We find that 
the flow speed peaks at |v| 1.0 in the reconnec¬ 

tion exhausts near the null, whereas farther out 
along the loop typical values drop to ~ 0.3. As 
the current layer rotates around the dome, sequen¬ 
tially reconnecting field lines, it produces a train 
of torsional Alfven waves that propagate along the 
loop (Fig. [^). Collectively, these waves form a 
large-scale traveling pulse whose angular velocity 
vector points along the loop toward the far-loop 
footpoints, until the leading waves arrive there 
and begin to reflect back toward the jet source. 
The reconnection proceeds relatively unhindered 
by these reflections, as the travel time to the far 
footpoint and back (ttravei) this case is roughly 
twice the duration of the jet (tjet)- Thus, we find 
that the returning flows arrive back at the dome 
well after the jet ceases. 

By comparison, the jet produced in the short 
loop {L/N = 1.46) starts much later {Brig ~ 840) 
and is much shorter in duration {tjet ~ 100). The 
later onset can be attributed to a combination of 
two factors: a significant fraction of the injected 
shear has already been transferred onto loop field 
lines via reconnection in the energy storage phase, 
and more twist is required to initiate the instabil¬ 
ity that drives the jet. The shorter duration of 
the jet is due primarily to the shorter travel time 
along the loop, which causes the reflections from 
the far loop footpoint to affect the jet source. We 
estimate the travel time along the loop and back 
to be tfravei ~ 50. The interaction of the jet re¬ 
gion and the return flows can be seen in Fig. [^-f. 
Following the onset of rapid energy release, fast 
reconnection-driven plasma flows are launched se- 


2010 2015 


10 














Fig. 6.— The impulsive jet produced wheu L/N = 2.40 (top) aud 1.46 (bottom) with maguetic field hues as 
iu Figure]^ Isosurfaces of velocity maguitude |v| = 0.14 are color-shaded accordiug to Vz, the out-of-plaue 
velocity compoueut, whose scale is saturated for clarity. The horizoutal-plaue shadiug shows as iu Figure 
Auimatious of both coufiguratious are available ouliue. 


queutially toward the apex of the loop, with the 
leadiug outflows traveliug aloug the shortest field 
hues uearest the photosphere (Fig. [^). Oue full 
travel time after ouset, at t = 900, the previously 
lauuched flows aloug the shortest field hues have 
reflected off the far loop footpoiut aud are returu- 
iug aloug the loop toward the recouuectiou regiou. 
Iu the iuterim, further outflows that curve over the 
top of this returu flow have beeu lauuched (Fig. 

Soou after, at t = 920, the jet outflows be¬ 
come fragmeuted (Fig.|^) as the couuterstreamiug 
flows iuteract. Both fast recouuectiou aud siguif- 
icaut maguetic euergy release theu cease as the 
recouuectiou at the dome is choked off by the re- 
turuiug flows. 

3.3. Relaxation Phase 

After the maiu euergy release phase has cou- 
cluded, both loops coufiuiug the jets relax toward 
a uew quasi-steady coufiguratiou. For L/N = 
2.40, the loug loop exteuds higher iuto the coroua 
where the field streugth is weaker aud the loops 
expaud more. As a result, the propagatiug twist 
compoueut of maguetic field expauds as it reaches 
the apex of the loop, theu uarrows agaiu as the 
disturbauce reaches the coujugate footpoiut ou the 
photosphere. The fastest flows occur iu a curtaiu- 


like baud arouud the periphery of the propagat¬ 
iug regiou of twist, where the sliugshot effect from 
the release of maguetic teusiou is strougest (Fig. 
&)• Couuterstreamiug flows aloug the loop are es¬ 
tablished ouce the leadiug wave reflects off of the 
far loop footpoiut. Iu additiou to the torsioual 
waves, a weaker lougitudiual oscillatiou of the loop 
is geuerated as the jet propagates aloug it, so that 
the loop sways relative to its liue-tied footpoiuts. 
The torsioual waves quickly distribute the iujected 
twist more eveuly aloug the loop as the separatrix 
dome relaxes (see auimatiou of Fig. [^-c), leav- 
iug the twisted loop to oscillate geutly as the as¬ 
sociated flows gradually damp away. Toward the 
eud of the simulatiou, weak curreut layers exteud 
aloug the expauded coroual loop (Fig. [^). These 
layers reflect mismatches iu ueighboriug flux tube 
leugths arisiug from the three-dimeusioual, iuho- 
mogeueous uature of the jet geueratiou. This pro¬ 
cess has beeu called recouuectiou-driveu curreut 
filameutatiou (Karpeu et al. 1996). A localized 
curreut layer also remaius at the uull, aud it cou- 
tiuues to slowly release the remuauts of twist from 
beueath the separatrix dome. 


For L/N = 1.46, the short loop with its re¬ 
duced expausiou iu height supports similar wave 
behavior, but iucreased dissipatiou from the couu¬ 
terstreamiug flows damps the waves more rapidly 
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Fig. 7.— The final state {t = 1200) in the two 
experiments, with magnetic field lines as in Figure 
The range and number of |J| contours has been 
increased to enhance the lower-amplitude features. 


(see animation of Fig. |^-f). In a manner simi¬ 
lar to the long loop, by the end of the simulation 
a large portion of the injected twist has been re¬ 
distributed along the loop, and extended current 
layers permeate the part of the loop affected by 
the jet. In this case, the currents in these layers 
are stronger, due to the shorter loop length over 
which the shear is spread (Fig. [^). A current 
layer also resides at the null in the final state. 

3.4. Reconnection Analysis 

To understand more quantitatively how recon¬ 
nection across the separatrix correlates with the 
observed jetting behavior, we investigated the con¬ 
nectivity of the magnetic field in the vicinity of the 
separatrix dome. Using our newly developed field- 
line integrating routine, we traced field lines from 
a 500^ grid of starting positions centered on the 
parasitic polarity. Each field line either starts and 


ends beneath the separatrix surface, or starts out¬ 
side of it and connects to a distant footpoint of 
the coronal loop. By labeling each starting point 
accordingly, one can see the footprint of the sep¬ 
aratrix surface at a given time. If reconnection 
occurs across the separatrix over time, this foot¬ 
print evolves. By tracking the evolution, the as¬ 
sociated reconnection rate may also be calculated 
(see Appendix [A|). 

Figure shows the footprint of the closed flux 
beneath the dome in each simulation, at four 
times: initially (a,e); at the onset of fast reconnec¬ 
tion (b,f); at the peak of the reconnection process 
(c,g); and in the post-jet relaxation phase (d,h). 
Figure shows the corresponding calculated re¬ 
connection rates. The rapid increase in reconnec¬ 
tion rate corresponds closely in each case to the 
onset of the jet in the volume, whereas the peak 
in reconnection rate occurs about halfway through 
the jet ( p.2[ ). 

For L/N = 2.40, it is clear that the separatrix 
surface has shifted only slightly during the energy 
storage phase (Fig. |^-b), consistent with the low 
reconnection rate during this period. At the on¬ 
set of the kink instability, the reconnection rate 
increases strongly (Fig. as the separatrix dome 
contorts. A channel of locally open field forms, 
penetrating deep into the previously closed field 
region as far as the polarity inversion line (Fig. 
[§). Once this channel is formed, the footprint 
of the surface rotates as highly sheared flux near 
the PIL is reconnected out of the dome and un¬ 
sheared loop flux is reconnected into it (see the 
online animation of Fig. [^. Because the field 
component normal to the photosphere in these 
simulations is held constant, the flux beneath the 
separatrix dome is a fixed quantity. To preserve 
it, the amount of flux being opened and closed 
across the separatrix must be equal at any given 
time (see also Pontin et al. 2013). This type of 
reconnection is often termed “interchange recon¬ 
nection.” The interchange of a section of highly 
sheared field with unsheared field reduces the av¬ 
erage shear of the closed flux beneath the dome, 
but does not remove it entirely. The reconnection 
in the later stages of the jet evolves the field to¬ 
wards a lower energy state that is nearer to poten¬ 
tial, and brings the majority of the recently opened 
(originally closed) flux back beneath the dome. In 
the final state, the separatrix remains somewhat 
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Fig. 8.— Footprint of closed field within the separatrix surface initially (a,e), at the time of fast reconnection 
onset (b,f), at the time of peak reconnection (c,g), and during the post-jet relaxation (d,f). White regions 
show magnetic flux that closes beneath the separatrix dome, whereas black regions show flux that connects 
to the distant coronal-loop footpoints. Red circles show the polarity inversion line in each case. Top panels: 
L/N = 2.40; bottom panels: L/N = 1.46. Y and Z are coordinate axes centered on the parasitic polarity. 
Animations of both configurations are available online. 



distorted because the closed field still retains part 
of the injected twist, whilst the remainder of the 
twist has been transferred to the enclosing coronal 
loop. At this stage, reconnection across the sepa¬ 
ratrix becomes very weak. Overall, the total flux 
beneath the dome is reconnected roughly twice - 
first when the sheared flux is opened, and second 
when it is closed again - with the opening and clos¬ 
ing occurring principally during the energy-release 
phase. 

By contrast, early reconnection in the more 
asymmetric configuration for L/N = 1.46 opens 
a channel into the sheared field next to the PIL 
during the energy-storage phase (Fig.|^). The re¬ 
connection rate ramps up at a nearly steady pace 
during the lead-up to the jet trigger time (Fig.[^. 
When the jet is triggered, the separatrix is already 
highly distorted and roughly all of the flux beneath 
the dome has been reconnected once ( ^4.2| ). The 
fast reconnection at the time of the jet (note the 
spike in Figure]^ rapidly closes down the previ¬ 
ously opened flux, returning it beneath the sepa¬ 


ratrix dome. Therefore, in this configuration the 
“opening” of flux occurs much more slowly during 
the energy-storage phase, whilst only the rapid re¬ 
closing of this flux occurs during the energy-release 
phase. This feature, in addition to the jet reflec¬ 
tion along the loop, helps to explain why this jet is 
so much shorter in duration than that in the long 
loop. 

4. Results II: Parametric Survey 

The results just described demonstrate that a 
domed 3D null point topology can produce jets 
in closed-field configurations. However, the rel¬ 
ative size of the jet source region {N) with re¬ 
spect to the enclosing coronal loop (L) plays an 
important role in determining the qualitative and 
quantitative features of the jet. Using the insights 
gained from the simulations shown above, we now 
describe the results of the entire parameter study. 


13 







1.0 



Fig. 9.— Reconnection rates calculated from the 
flux swept out by the separatrix surface. The 
times of the connectivity maps shown in Figure 
are marked by asterisks (reconnection onset; Fig. 

3,f) and diamonds (peak reconnection rate; Fig. 


4.1. Initiation 

We begin by discussing the conditions for jet 
onset in the different configurations. All but one 
of the cases that we studied eventually produced 
an impulsive jet, although some were very weak. 
We discuss the exception (with the largest value 
of L/N) at the end of this section. All jets oc¬ 
curred during either the ramp-down phase of the 
driving or after the driving ceased. We defined the 
time at which each jet was initiated {ttrig) as the 
time when the rate of magnetic-energy liberation 
suddenly increases (details are discussed below in 


Figure 10 r shows how ttrig varies with L/N. 
For the largest values of L/N, the null point is 
nearest the top of the dome and the dome as a 
whole has the greatest cylindrical symmetry (cf. 
Fig. [^). This symmetry inhibits the initiation 


of the kink instability that drives the jet (Pariat 
et al.||2QQ9[ |Rachmeler et ^|2Q1Q ). Consequently, 
as L/N increases the jet is increasingly delayed 
and ttrig increases. For the smallest values of L/N, 
the null point is farthest to the side of the dome 
and the dome as a whole has the greatest asym¬ 
metry (cf. Fig. 1^). Current-sheet formation and 
reconnection occur more readily at the null, reduc¬ 
ing the rate of energy buildup beneath the dome 
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Fig. 10.— For each simulation vs. aspect ratio 
L/N: (a) Jet trigger time, ttrig] (b) At t = Urig] 
injected relative helicity normalized by the square 
of the flux beneath the dome, Htrig (diamonds), 
and the measured maximum number of turns, 
Mtrig (asterisks). 


(^3.4). This effect delays the jet, so that ttrig also 
increases as L/A^ decreases. The shortest trigger 
time occurs at an intermediate value of the aspect 
ratio, L/N 2.1, where the two delaying effects 
are jointly minimized. 

Similar trends are seen in the maximum number 
of turns, M, present in the closed field at t = ttrig 
{Mtrig = M{ttrig)) aud iu the normalized relative 
helicity, injected into the volume by this time 
{Htrig = Hn{ttrig)), Fig.[^. The former is calcu¬ 
lated by counting the number of turns about the 
center of the parasitic polarity for each of the 500^ 
field lines traced to analyze the magnetic connec¬ 
tivity. The latter quantity is calculated from 


Hn{t) = ^ 


(t) 


dome 


( 8 ) 
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where the injected helicity is 

m = f 
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dt' 


dH, 


mj _ 

dt 


j (v • A) B:j,dS (9) 

2 j {vyAy B^dydz. 

The last integral, evaluated over the photosphere 
{x = 0), is derived from the |Finn & Antonsen 


( 1985| ) gauge-invariant form of the helicity, 
H{t) = y (A + Ap) ■ (B - Bp) dV. 


( 10 ) 


Here Ap and are the current-free magnetic 
potential and field, respectively, for the instanta¬ 
neous Bx distribution at the photosphere. The 
normalized helicity, measures the average 

twist injected into the flux beneath the separa- 
trix dome. Its value when the jet is initiated, 
Htrig, can be seen to match very well the max¬ 
imum number of turns, Mtrig^ measured in the 
volume numerically at this time. The two devi¬ 
ate for the shortest, most asymmetric systems, 
where reconnection occurs across the initial sep- 
aratrix dome prior to jet onset. For L/N > 2.1, 
the increasing cylindrical symmetry allows more 
twist to be stored, increasing the number of turns 
and the helicity injected into the field before the 
jet is triggered. For L/N < 2.1, the delayed jet 
trigger allows more helicity to be injected and 
more twist to be imparted to the field. However, 
the increasingly important spine-fan reconnection 
in this lower range spreads some of the injected 
helicity along the coronal loop and restricts the 
strongly twisted field lines to a smaller portion of 
the flux beneath the dome ( p.4[ ). The helicity 
eventually plateaus at small T/A, whilst the num¬ 
ber of turns peaks and then declines again as the 
reconnection penetrates farther toward the PIL. 
The turn values that we measured lie within the 
range Mtrig ^ [0.7,1.1], somewhat below the 1.4 
turns required for the initiation of the kink insta¬ 
bility in open-field configurations with perfect ini¬ 
tial symmetry ( Pariat et al.|2QQ9 Rachmeler et al 
2010) and consistent with the 0.8 turns required 


in a setup with a tilted background field (Pariat 
et al.|201Q ). This supports our interpretation that 
the impulsive energy-release phase is driven by the 
onset of a kink-like instability. 


Similar results for the variation in jet trigger 
times for magnetically open configurations have 
been presented by Pariat et al. (2015). They in¬ 


vestigated the influence of the inclination angle of 
a straight, uniform, open background field on the 
initiation and evolution of jets. Their simulations 
showed the same two effects governing the trig¬ 
ger time of open jets: the kink instability of the 
twisted field, which occurs later in configurations 
with more vertical background fields and greater 
cylindrical symmetry; and spine-fan reconnection 
across the separatrix, which lowers the rate of en¬ 
ergy storage beneath the dome and delays the jet 
in configurations with more horizontal background 
fields and greater dome asymmetry. The shortest 
trigger time occurred for an intermediate inclina¬ 
tion angle of the field. 

However, there are two important differences 
between those open-field investigations and our 
closed-field simulations. First, our background 
magnetic field falls off with height above the pho¬ 
tosphere. Consequently, the separatrix domes in 
our simulations expand essentially vertically (cf. 
Figs. [^,f), rather than along the direction of the 
uniform background field as occurs in the open 
cases. Second, as the dome expands into the closed 
loop, the tilt angle of the background field at the 
null increases as the field drapes over the dome. 
This generates a strengthening non-axisymmetric 
current layer at the null in all configurations, lead¬ 
ing to initially weak reconnection outflows that in¬ 
crease as the dome expands. Therefore, there is no 
completely axisymmetric configuration as can be 
achieved in the open-field case, where reconnec¬ 
tion is geometrically inhibited during the buildup 
phase. Given these differences, it is not surpris¬ 
ing that the trigger-time-minimizing angle in our 
simulations differs from the ~ 8° tilt angle identi¬ 


fied by Pariat et al. (2015). Although the effective 
inclination angle in our simulations changes with 
time, a rough estimate of the initial angle (^o) can 
be obtained by using the null position relative to 
the center of the parasitic polarity at t = 0. This 
is a close approximation to the angle of the spine 
lines near the null, and is equivalent to the angle 
of the straight background field in open configu¬ 
rations. An approximately inverse linear relation¬ 
ship between Oq and L/N is found, as shown in 
Figure [m The aspect ratio that minimizes the 
jet trigger time for our closed jets, L/N 2.1, 
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Fig. 11.— For each simulation vs. aspect ratio 
L/N: Inclination angle, Oq. 

corresponds to a relatively steep inclination angle, 
(9o ~ 20°. 

We conclude this subsection by discussing 
briefly our longest-loop configuration {L/N = 
2.73), which remained stable and did not produce 
an impulsive jet. This case is most cylindrically 
symmetric, with the null point remaining near the 
top of the dome throughout the simulation. The 
number of injected turns {Mt^ig ~ 1-2) was less 
than the critical amount {Mtrig ~ 1-4) needed 
to set off the kink instability in an open field 
with zero inclination angle. Draping of field lines, 
mentioned above, formed a thin current layer sur¬ 
rounding the single null point. The reconnection 
associated with this current layer was very weak, 
did not reach the sheared field beneath the dome, 
and did not destabilize the configuration. We ran 
this simulation 400 time units past the end of the 
driving period, during which time the magnetic 
field relaxed and reconnection within the current 
layer tapered off almost entirely. The final state 
contains a current layer around the null and, evi¬ 
dently, is stable. 

4.2. Energetics 

We now consider the energetics of our jets. Be¬ 
cause each configuration has a slightly different 
magnetic flux and peak field strength associated 
with the parasitic polarity, the energy injected by 
the boundary driving is different in each case. The 
total energy injected into the magnetic field by the 
footpoint motions up to time t is calculated by 
integrating the Poynting flux across the bottom 


^inj{t) — I dt' j c(E X B) • ndS 

= - J J J +VzBz) Bx dy dz. 

( 11 ) 

To compare the runs equitably, we normalize the 
various energies with respect to the total energy 
injected into the volume by the footpoint motions 
in each case, Ej/Zj = Einj{t oc). Table lists 
the resulting values of El//j for each configuration. 

Early in each simulation, essentially all of the 
injected energy is stored in the coronal magnetic 
field, as very little is converted to kinetic or ther¬ 
mal energy of the plasma. Because the component 
of the magnetic field normal to the photosphere is 
held fixed throughout each run, the current-free 
field Bp and its energy Ep also are independent 
of time. Thus, the instantaneous free magnetic 
energy is given simply by 

^Emag{t) ^mag {t)-Eraaa{Q). ( 12 ) 

Figure [T^ shows the evolution of the normal¬ 
ized AEmag{t) cach of our jet-producing simu¬ 
lations. Also shown is the profile of Einj{t), which 
after normalization is nearly identical for all cases. 
The onset of the jet in each configuration is sig¬ 
naled by a rapid decrease in AE^ag • As previously 
mentioned, jet onset occurs over a broad range of 
times with respect to the phase of the boundary 
driving: some occur during the ramp-down phase 
(500 < t < 1000), others after the driving ceases 
{t > 1000). To characterize the time and duration 
of each jet, we focused on the cumulative energy 
liberated from the magnetic field, 

Eiih{t) — Einj (t) - AEmag{t). (13) 

This quantity, shown in Figure \T2\ ) after being 
normalized to measures the total energy re¬ 

leased by the magnetic field whilst taking into ac¬ 
count the different absolute energy injection rates 
at the different times of jet onset. We define the 
trigger time, Lrig, as the time when a noticeable 
increase in Em^ occurs (corresponding to fast en¬ 
ergy release as the jet begins) and the jet dura¬ 
tion, tjet, as the time elapsed thereafter until the 
sharp increase in Em^ subsides. Table lists Lrig 
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Fig. 12.— Energies normalized to the total injected energy, vs. time t for each simulation: (a) Injected 
energy Einj{t) (black lines) and stored magnetic energy AEmag{t) (colored lines); (b) Liberated magnetic 
energy Eiii){t); (c) Kinetic energy Ekin{t)- In (a)-(c), diamonds and triangles mark the values at t = ttrig 
and t = ttrig + tjet^ respectively. Em, vs. aspect ratio L/N for each simulation: (d) Liberated energy up to 
time t = ttrig {Etrig^ boxcs) and liberated energy during the energy release phase 0 < t — ttrig < {Ejet^ 
asterisks). 


and tjet in each simulation. Diamonds and trian¬ 
gles in Figure [T^ mark the normalized Em, values 
at times ttrig and ttrig + tjet^ respectively. These 
times are well-correlated with the rapid increase 
and subsequent decline of the volumetric kinetic 
energy, Ej^in{t), following the launch, travel, and 
deceleration of the jet in each case (Fig. [T^). 

The transition in behavior around L/N = 2.1, 
from jets preceded by significant reconnection in 
the energy buildup phase to jets with relatively 
little, is evident in the magnetic energy curves. In 
Figure [T^, the injected energy closely matches the 
free magnetic energy prior to the jet {t < ttrig) for 
configurations with L/N > 2.1 (dashed lines), due 
to the weak reconnection occurring in this phase. 
Correspondingly, the liberated energy Em, is small 


until the onset of the impulsive phase of these jets, 
when large increases in Eut, are evident (Fig.[T^). 
In contrast, for L/N < 2.1, the injected and free 
energy curves deviate early in the evolution (Fig. 


12r, colored solid lines), with significant fractions 


of the injected energy liberated prior to jet on¬ 
set in the most asymmetric (smallest L/N) con¬ 
figurations (Fig. [T^). Despite the marked differ¬ 
ences in the timing of energy liberation among the 
various cases, a rather consistent total of 50% to 
60% of the injected energy is liberated during the 
combined energy-storage and -release phases (Fig. 


Hit). 

Figure [T^ shows that the energy liberated by 
time t — ttrig 1 


Etrig Em,{ttrig) 1 


(14) 


17 

























normalized to is as large as 20% for the cases 
with small ratios oi L/N and drops to less than 
5% for configurations with L/N > 2.1. The slight 
increase at the largest values of L/N is due to 
late-time, weak reconnection associated with drap¬ 
ing of field lines over the strongly expanded dome. 
The energy liberated during the impulsive energy- 
release phase itself, 

Ejet T 5 

normalized to increases from about 40% for 
moderate ratios of L/N to 60% for the largest. 
There is a sharp fall-off in this quantity at small 
aspect ratios {L/N < 1.3), due to the signifi¬ 
cant reconnection in the energy storage phase and 
the jet reflection along the short loops for these 
highly asymmetric cases. The impaired energy re¬ 
lease produces weak jets with small kinetic ener- 
gies (Fig.[l^). 

A deeper understanding of the energetics of our 
jets results from considering the cumulative recon¬ 
nected flux, normalized to the flux within the sep- 
aratrix dome, shown in Figure [T^. The resem¬ 
blance between the curves of free-energy libera¬ 
tion (Fig.jl^) and of interchange-reconnected flux 
across the separatrix (Fig. [T^) is striking. Close 
similarities can also be seen between the magnetic 
flux reconnected prior to each jet, 

^trig — ^rec{^trig^ •) 

and during the energy-release phase. 


other characteristic jet quantities, the transition 
between these two behaviors of the reconnected 
flux occurs at L/N ^ 2.1. 

These trends in the energy liberated and the 
flux reconnected also are reflected in the dura¬ 
tions of the jets and in the associated peak kinetic 
energies (Fig. [T^,b). Those configurations with 
large values of L/N^ which reconnected the most 
flux and liberated the most free energy during 
the energy-release phase, produced the longest- 
lived jets with the greatest kinetic energies (up 
to 10% of the total injected energy), whilst 
those with the smallest ratios oi L/N produced 
the shortest-lived, least-energetic jets. The dura¬ 
tion almost triples in length, from tjet ~ 120 to 
~ 340, with a steep transition around L/N = 2.1 
(Fig. pA^), correlated with the reduction in nor¬ 
malized flux reconnected prior to jet onset (Fig. 
[T^ ; boxes). The peak kinetic energies show a 
more complex dependence upon the aspect ra¬ 
tio L/N^ but generally increase for larger ratios. 
Jets in shorter loops (small L/N) have signifi¬ 
cantly shorter travel times along the loop {Uravei) 
than the duration of the jet itself {tjet)i some¬ 
times yielding double-peaked kinetic energy curves 
as the counterstreaming flows interfere with each 
other. For higher values oi L/N the travel time 
increases, but this is mediated in part by the in¬ 
crease in jet duration. By inspection, we And that 
jets with L/N > 1.7 are relatively unaffected by 
reflections from the far loop footpoints, although 
they do exhibit counter streaming flows in the loop. 


^jet ^recij^trig T ^jet) ^rec {ttrig). (17) 


5. Correspondence with Observed Jets 


normalized to ^dome^ shown in Figure [T^, with 
the corresponding liberated energies over these 
time intervals (Fig. [T^). Both show the direct 
link between the rates of interchange reconnection 
and energy release prior to and during the jet. 


By the end of the energy-release phase, all but 
the most compact, asymmetric configurations re¬ 
connect the flux beneath the dome at least twice 
(4^rec > 2Tc^ome)- oncc when the sheared held 
near the PIL reconnects to distribute twist along 
the coronal loop, and again when it reconnects 
back down beneath the dome. However, as dis- 
for small L/N the opening oc¬ 


cussed in ^3.4 


curs during the energy-storage phase, whilst for 
large L/N the opening and closing are delayed 
until the energy-release phase. As was found for 


5.1. Qualitative Features 

In all of the configurations that we studied, the 
jet produced was confined by a coronal loop whose 
footprint envelopes the jet source region. This is 
a generic feature of all closed-field jets, most evi¬ 
dent when the coronal loop has significant curva- 


ture (e.g. Shibata et al.|1992[ Shimojo et al. 

1998 

Torok et al.|2009| Yang et al.|2012 Guo et al. 

2013 

Lee et al.|2013t Schmieder et al.|2013t Zheng 

et al. 

2013| Cheung et al.||2015). The curtain-like shape 


of our jets was shown to expand as it propagates 
towards the loop top (driven by the expansion of 
the ambient held as the background held strength 
declines with height) and contract again toward 
the far loop footpoint. This effect is expected to 
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Fig. 13.— Cumulative reconnected fluxes, ^rec, normalized to the flux beneath the separatrix dome, 
for each simulation: (a) Evolution vs. time t, with values marked when t = ttrig (diamonds) and t = ttrig-\-tjet 
(triangles); (b) Variation vs. aspect ratio L/N at time t = Lrig {^trig^ boxes) and during the energy-release 
phase 0 < t — Lrig < jeti astcrisks). 
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be most observable in large jets that have a rela¬ 
tively wide source region, so that the affected loop 
strands expand noticeably in the corona and the 
jet flows are energetic enough to defy gravity and 
reach high up along the loop. The animation of 
the active region jet in Fig.l^shows such a large jet 
where this expansion and contraction are clearly 
visible. 

Our model predicts differing behaviors for the 
waves launched during the impulsive jet phase de¬ 
pending upon the duration of the jet (tjet) com¬ 
pared with the travel time along the coronal loop 
i/^travei)' When L/N is very large, ^jet 

and the jet-launched wave motions are expected to 


travel freely along the coronal loops. [Torok et al. 


(2009) observed and modeled an example of such 
freely propagating torsional wave motions along 
longer coronal loops. When the two time scales 
are comparable, our model predicts that the jet 
generation is relatively unhindered by reflections 
from the far loop footpoint, but counterstreaming 
flows will be present in the loop. A possible mani¬ 
festation of this was described by Qiu et al. (1999), 
who observed counterstreaming flows in Ha loops 
in an active region at the same time that flaring 
and jetting were observed in a mixed-polarity re¬ 
gion at one end of the loop system. When L/N is 
small and Uravei < our model predicts that 
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the reflected jet-launched waves will interfere with 
the further development of the jet. The interac¬ 
tion of the returning flows and the jet outflows 
could conceivably generate turbulence in the loop, 
leading to extended emission along the connecting 
loops as well as within the anemone region. |Shi- 


mojo et al. (1998) described soft X-ray emission 


spread along short coronal loops where the esti¬ 
mated travel time was shorter than the jet lifetime, 
lending some credence to this idea. 

As in open-field regions, closed-field jets some¬ 
times are generated repeatedly from the same 
source region (e.g. Cheung et al.|[2M5 ). The free 
energy that drives the jets in our model arose from 
boundary motions, applied on the photosphere, 
which were stopped once a single jet was pro¬ 


duced. In open fields, Pariat et al. (2010) showed 


that homologous jets are generated by maintaining 
the driving. We tested this on one of our simula¬ 
tions by maintaining the photospheric driving, and 
found that it too produced homologous jets along 
the same coronal loop. Thus, our model can ex¬ 
plain the observed homology of some closed-field 
jet regions. 

5.2. Quantitative Measures 

Quantitatively, the properties of our jets can 
be compared to those observed by applying scale 
factors ps, and Lg to the dimensionless sim¬ 
ulated quantities to obtain typical coronal val¬ 
ues of mass density, magnetic field strength, and 
length, respectively. The associated scale factors 
for pressure (Ps), velocity (Vg), time (t^), energy 
(Eg), and magnetic reconnection rate (flux per 
unit time; Tg) are 

Ps = Bl 

v=E. 

ts = (18) 

Es = BlLl 

ifg = BgVgLg. 

For simplicity, in this discussion we will work 
solely with orders of magnitude in the scale fac¬ 
tors. 

The most elementary example is jets occurring 
in areas of quiet Sun {qs). There, we can assume 


Bg = hence the strengths of the coronal-loop 
{Bel) and parasitic-polarity (Bpp) fields are 


Bit 


4G, 


Bi^ 

pp 


21 G. 


(19) 


The pressure scale factor = 1, so the thermal 
pressure {Pth) is 


^th 


1 X 10 ^ dyn cm 


( 20 ) 


At a temperature of 1 x 10^ K, the corresponding 
mass density is 


^ 1 X 10“^^ g cm 


-3 


( 21 ) 


Hence, in quiet Sun we set p^ = 1 x 10“^^, which 
together with Bg = 1 gives Vg = lx 10^. Our peak 
dimensionless jet velocity is ~ 1.0, so we find 


vf:, 


lxl0®cms“r (22) 

This is consistent with the maximum apparent jet 


flow speeds reported by Shimojo et al. (1996) and 


Savcheva et al. (2007). Typical velocities within 


our main curtain-like spray were somewhat lower, 
with dimensionless values ~ 0.3, scaling to 3 x 10^ 
cm s“^. This falls within the ranges of velocities 


reported by Shimojo et al. (1996) and Savcheva 
et al. (2007), and is consistent with the jet speeds 


reported by Shibata et al. (1992). Finally, a length 
scale factor Pg = 1 x 10^ applied to our average 
dimensionless loop length gives a coronal value 


^ 12P, 1.2 X 10^^ cm. 


(23) 


For this long coronal loop, the remaining scale 
factors are p = 1 x 10^, Eg = 1 x 10^^, and 
Tg = 1 X 10^^. Using values from Table [ij for 
the typical jet duration, peak kinetic energy, and 
peak reconnection rate we obtain, respectively. 


pet " 

s 200G 



s 2 X 10^ s, 

(24) 


s .04 X 125Es 



a 5 X 10^^ erg, 

(25) 

^ rec 

a 1.04', 



a 1.0 X 10^^ Mx s“\ 

(26) 


The duration of our quiet-Sun jet is roughly the 
mean of the large ran ge of jet lifetimes reported by 
Shimojo et al. (1996) (?^ 2 x 10^-2 x 10^ s) and is 
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at the high end of the range reported by [Savche^ 
et al. (2007) 2 x 10^-2 x 10^ s). Its kinetic en- 


erg y is near the upper e nd of the range reported 
by Shibata et al. (1992) 1 x 10^^ erg). The 

reconnection rate varies from about 1 x 10^^ Mx 
s“^ during the early, slow-reconnection (energy- 
storage) phase to 1 X 10^^ Mx s“^ during the im¬ 
pulsive, fast-reconnection (energy-release) phase 
of the jets. 

For jets occurring in active regions (ar), we set 
Bs = 10. The strengths of the coronal-loop {Bd) 
and parasitic-polarity (Bpp) fields then are 


B^[ ^ 40 G, B^^ ^ 210 G. (27) 

The pressure scale factor Pg = 100, so now the 
thermal pressure is 

« 1 X 10° dyn cin“2. (28) 

Here, the corresponding mass density is 

1 X 10“^"^ g cm-^ (29) 


We therefore set pg = 1 x 10“^^, which together 
with Bg = 10 again gives Vg = 1 x 10^. Thus, as 
in quiet Sun we find 

1/“^ « 1 X 10® cm s“r (30) 

Here, we assume a length scale factor = 1 x 10® 
to model jets in compact active-region loops, 

L“'' « 12L, ^ 1.2 X 10° cm. (31) 


The rest of the scale factors are tg = 1, Eg = 
1 X 10^^, and (as before) = 1 x 10^^. For the 
typical jet duration, peak kinetic energy, and peak 
reconnection rate we obtain, respectively. 


Get " 

s 200ts 



s 2 X 10^ s, 

(32) 


s .04 X 125E, 



a 5 X 10^° erg, 

(33) 

rec 

a 1.04-, 



a 1.0 X 10^^ Mx s“\ 

(34) 


The duration of our active-region jet is at the low 
end of the range of jet lifetimes reported by |Shi - 


mojo et al. (1996) and Savcheva et al. (2007), and 


somewhat shorter than the observed jet shown in 


Figure Its kine tic energy is near the mean of the 
range reported by Shibata et al. (1992) (^ 1 x 10^^ 
- 1 X 10^^ erg). The reconnection rate varies over 
the same range as in quiet Sun, 1 x 10^° Mx s-i 
to 1 X 10^^ Mx s“^. 


6. Discussion 


We have investigated the initiation and evo¬ 
lution of solar jets in closed coronal loops via 
three-dimensional numerical simulations of the 


embedded-bipole model (Antiochos 1996), which 
has been investigated extensively to explain 


coronal-hole jets (Pariat et al. 2009 2010 2015) 


As in those prior studies, our impulsive jets con¬ 
sist of twisted, curtain-like outflows that exhibit 
strong helical motions and filamentary structure. 


resembling observations (e.g. Torok et al. 

2009 

Yang et al.||2012t 

Guo et al. 1120131 Lee et al. 

2013 

Schmieder et al.| 20131 |Zheng et al. 2013| Ghe-| 

ung et al. 2015). 

Due to the variation in loop 


cross-section with distance along the loop, our jets 
expand and then contract laterally as they prop¬ 
agate along the loop toward its far footpoint, as 
observed. By scaling our dimensionless results us¬ 
ing typical input parameters characterizing coro¬ 
nal jet sources ( Shibata et al.|19^ Shimojo et al 
1996), we obtain values for the jet speeds, dura¬ 


tions, and energies that are consistent with those 
observed (^5.2|). 


Our results further show that the evolution is 
highly sensitive to the relative sizes of the closed 
dome (N) of the jet source region and the coronal 
loop (L) within which the source is embedded. We 
found that configurations with large L/N ratios 
store the greatest amount of magnetic free energy 
and produce the longest-duration, most energetic 
jets; those with small L/N ratios release a signif¬ 
icant fraction of the injected free energy prior to 
onset of their shorter-duration, less energetic jets. 
The transition between these behaviors occurs at 
L/N ^ 2.1, where the number of turns of induced 
twist required to initiate the impulsive jet is min¬ 
imized (Fig. [T^). Smaller configurations liberate 
up to 20% of the stored free energy prior to jet 
onset; larger ones liberate essentially none (Fig. 
[T^). These energies reflect amounts of cumula¬ 
tive reconnected flux prior to onset of as much as 
100% and as little as 10%, respectively, of the total 
flux enclosed beneath the dome (Fig. [T^). Dur- 
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ing the fast reconnection that drives the impulsive 
jet, in contrast, the liberated energy ranges from 
less than 10% for small L/N to as much as 60% 
for large L/N (Fig.[T^), and the reconnected flux 
ranges from as little as 30% to as much as 200% 
(Fig.[^). The corresponding jet durations range 
from fewer than 100 to more than 300 Alfven times 
(Fig. [T^). Finally, the kinetic energies in the 
jet flow, normalized to the total injected energy, 
range smoothly from essentially zero at L/N ^ 1 
to about 5% at L/N ^ 2.4, and then rise at a 
more rapid pace toward higher L/N (Fig.pTjp). In 
configurations with L/N > 2.4, the reconnection 
driving the impulsive jet ceases before the jet flow 
reaches the far footpoint of the loop and reflects 
back into the jet source region. In configurations 
with L/N < 1.7, on the other hand, the coun¬ 
terstreaming forward and backward (reflected) jet 
flows interfere strongly with each other (Fig. [^, 
and even appear to choke off the reconnection out¬ 
flow in our most compact configurations. 

Following the cessation of the impulsive recon¬ 
nection and subsidence of the principal jet out¬ 
flows, our system relaxes toward a new quasi¬ 
steady state with filamentary current structures 
threading the coronal loop. These structures arise 
from a local mismatch of neighboring flux-tube 
lengths, driven by the spatially and temporally 
intermittent transfer of twist to the loop during 
the three-dimensional evolution: reconnection- 
driven current filamentation ( Karpen et al.||1996 ). 
At higher grid resolution, these currents should 
become even more filamentary in structure and 
greater in strength, and could produce quasi¬ 
steady heating of the loop in the aftermath 
of the jet. In addition, it is known that 3D 
null-point current layers are explosively unsta¬ 
ble to resistive tearing at high Lundquist numbers 
(Wyper & Pontin 2Q14a|b ). Attaining the requi¬ 
site Lundquist-number threshold demands better 
resolution than the simulations undertaken in this 
work. Sufficiently high-resolution simulations also 
should reveal more fine structure in the jet itself, 
plausibly including the formation and ejection of 
small-scale plasma concentrations (“blobs”; Zhang 


& Ji 2014) in the jet outflow. These intriguing pos¬ 


sibilities are being evaluated in our ongoing study 
of coronal loop jets. 
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A. Reconnected Flux and Reconnection Rate 

To calculate the reconnection rate in our simulations we take advantage of the fact that there is a true 
separatrix surface across which the rate of flux transfer may be measured. In each of the simulations the 
boundary driving is localized within the circular polarity inversion line so that near the separatrix the field 
line footpoints are line tied and fixed in position. Any change in the position of the separatrix is therefore 
due to reconnection occurring within the volume. Since the driven ends of the closed field lines always 
remain beneath the dome we can ignore the fact that these field lines are moving relative to the position of 
the starting grid from which the field lines are traced. 

We calculate the instantaneous reconnection rate in the following manner. For each field line in our tracing 
grid we assigned a magnetic flux element 


Ay Adz, (Al) 

where Ay and Az are the separation of the starting positions in the grid and {Bn)ij is the magnetic field 
normal to the photosphere at this position. At each time we compare the connectivity of each field line with 
the previous time. The flux elements of all the field lines which have been “opened” and “closed” in this 
time interval are summed to give the total opened and closed flux respectively 


opened — E elosed — E 

opened i,j elosed i,j 


(A2) 


To a high degree of accuracy in each of our simulations, the opened flux matched the closed flux. The 
reconnection rate is taken to be the simple average 


^( 4 ) 


opened T elosed 

2A4 


(A3) 


Here tp is the average time of evaluation and Atk is the time increment, 
time t = is the sum 

K 

^ree (t) = E^(4)A4. 

k=l 


The flux reconnected up to any 


(A4) 
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